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We describe a method for comparing the frequency spectra of oscillating stars. We focus on solar-like 
oscillations, in which mode frequencies generally follow a regular pattern. On the basis that oscillation 
frequencies of similar stars scale homologously, we show how to display two stars on a single echelle 
diagram. The result can be used to infer the ratio of their mean densities very precisely, without 
reference to theoretical models. In addition, data from the star with the better signal-to-noise ratio can 
be used to confirm weaker modes and reject sidelobes in data from the second star. Finally, we show 
that scaled echelle diagrams provide a solution to the problem of ridge identification in F-type stars, 
such as those observed by the CoRoT space mission. 

1 Introduction 

This paper discusses how to compare frequency spectra of oscillating stars. We focus on solar-like 
p-mode oscillations, in which mode frequencies generally follow a regular pattern. This makes it useful 
to characterize them by a handful of frequency separations: the so-called large separation Ais between 
consecutive overtones of a given angular degree I, and the small separations between adjacent modes 
of different degree. These frequency separations have the advantage of being closely related to physical 
properties of the stellar interior (see Section 2). Measuring them and their variations with frequency 
and comparing with theoretical models is a major focus of asteroseismology (see, for example, reviews 
by Brown & Gilliland 1994 and Christensen-Dalsgaard 2004). 

There is somewhat less focus on absolute frequencies, partly because stellar models do not properly 
model the near-surface layers (Christensen-Dalsgaard et al. 1988; Dziembowski et al. 1988; Rosenthal 
et al. 1999; Li et al. 2002). This makes it difficult to compare individual observed frequencies with 
models, although Kjeldsen et al. (2008) have proposed an empirical correction that appears promising, 
at least for stars reasonably close in effective temperature to the Sun. 

Here, we wish to compare observations of one star with observations of another, and so difficulties 
with models are not relevant. We are motivated by the expectation from homology that if two stars 
are sufficiently similar, their oscillation frequencies will be in the same ratio as the square roots of their 
mean densities: 



Here, we are comparing modes in the two stars with the same radial order (n) and angular degree (I). 
Even if the two stars are not similar, we might still expect Equation 1 to provide a useful approximation. 
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Of course, it also follows that the large separation scales in the same way: 



Av 2 



Ai/i 




(2) 



We now present some examples and applications, using echelle diagrams to visualize the comparisons 
between stars. 



The echelle diagram, first introduced by Grec et al. (1983) for global helioseismology, is nowadays 
used extensively in asteroseismology as a valuable way of displaying oscillation frequencies. It involves 
dividing the spectrum into segments of length Av and stacking them one above the other so that modes 
with a given degree align vertically in ridges. Any departures from regularity, such as variations in the 
large separation with frequency, are clearly visible as curvature in the echelle diagram, and variations in 
the small separations appear as a convergence or divergence of the corresponding ridges. 

We conventionally define three observable small frequency separations: Sv 2 is the spacing between 
/ = and I = 2; <5i^i3 is the spacing between I = 1 and I = 3; and Svoi is the amount by which I = 1 is 
offset from the midpoint of the I — modes on either side. In practice, the large and small separations 
are observed to vary with frequency. 

The regularity in solar-like oscillation spectra allows us to write the mode frequencies in terms of 
the large and small separations, as follows: 



where e is a dimensionless offset. The small separation di is zero for I = (radial modes), and equals 
5voi for I = 1, Sv a2 for I = 2 and (Suoi + Svi 3 ) for I = 3. Bedding et al. (2010b) have suggested that 
this last separation should be designated Sv ns . 

Equation 3 describes the oscillation frequencies from an observational perspective. A theoretical 
asymptotic expression (Tassoul 1980; Gough 1986 2003) relates Av, di and e to integrals of the sound 
speed. In particular, Av measures quite accurately the sound travel time across the diameter of the 
star, while the small separations are sensitive to the structure of the core and e is sensitive to the surface 
layers. 

When making an echelle diagram, it is common to plot v versus (v mod A^), in which case each 
order slopes upwards slightly. However, for grayscale images it can be preferable to keep the orders 
horizontal. We have done that in this paper, and so the value given on the vertical axis is actually the 
frequency at the middle of the order. 

3 Scaled echelle diagrams and their applications 

The ridges in an echelle diagram will only appear vertical if we use the correct value of Av. For this 
reason, it does not generally make sense to plot two stars on the same echelle diagram. However, if 
the frequencies of the second star have been scaled by multiplying them all by the ratio of the large 
separations, we are led by Equation 1 to expect that its ridges can be made to coincide with those of 
the first. The scaling factor can be fine-tuned to optimize the alignment in two different ways: 

1. by matching the slopes for the two stars (making both vertical), which means matching Av, or 



2 Echelle diagrams and the asymptotic relation 



v n j = Av(n + 



\l + e) - di 



(3) 



2. by overlaying the ridges as closely as possible, although this may mean they have different slopes, 
which means forcing them to have the same value of e. 
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Figure 1: Scaled echelle diagram comparing three main-sequence stars. The greyscale is the power spectrum of 
q Cen A, filled symbols are oscillation frequencies for the Sun after multiplying by 0.7816, and open symbols are 
frequencies for a Cen B after multiplying by 0.6555 (see text for references). The scaling factors were fine-tuned 
to make the ridges for all three stars parallel (Method 1). Symbol shapes indicate mode degree: I — (circles), 
I = 1 (triangles), I — 2 (squares) and I = 3 (diamonds). 



Method 1 is shown in Figure 1 for three main-sequence stars: a Cen A (Av ~ 106 /xHz), the Sun 
{Av w 135 £iHz) and a Cen B (Av « 162 /U.Hz). The greyscale is the power spectrum of a Cen A 
as observed in a two-site campaign with UVES and UCLES (Bedding et al. 2004), but with weights 
optimized to minimize the sidelobes (Arentoft et al. 2010). The filled symbols are oscillation frequencies 
for the Sun (Broomhall et al. 2009; Table 2) after multiplying by 0.7816, and open symbols are 
frequencies for a Cen B (Kjeldsen et al. 2005) after multiplying by 0.6555. The scaling factors were 
tuned (using simple 'trial and error') to make the ridges for all three stars parallel. Note that scatter 
of observations about the smooth ridges for a Cen A and B is due to the relatively short duration of 
the observations (only a few times longer than the mode lifetimes). 

There is a systematic progression in stellar parameters (mass, effective temperature and luminosity) 
as we go from a Cen A through the Sun to a Cen B. We see in Fig. 1 a corresponding progression in 
the positions of the ridges, which corresponds to a change in e (see Equation 3). Apart from this, there 
is a close similarity between the three stars, although there are subtle differences in the curvatures of 
the ridges and in the small separations between them. 

If two stars have very similar parameters, the ridges will almost coincide (Methods 1 and 2 become 
the same). This is the case for the Sun and the solar twin 18 Sco, recently observed with HARPS and 
SOPHIE (M. Bazot et al., in prep.), for which the scaling factor gives an extremely precise measurement 
of the mean density. In this case, we can also use one star as a guide when identifying modes in the 
other (and eliminating aliases). Another example is shown in Figure 2, for a pair of low-mass stars: the 
filled symbols show observed frequencies for r Cet (Av « 170 ^Hz; Teixeira et al. 2009) and the open 
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Figure 2: Scaled echelle diagram comparing two simliar stars low-mass stars. Filled symbols: frequencies 
(unsealed) for r Cet. Open symbols: frequencies for a Cen B multiplied by 1.0478. The scaling factor was 
tuned to make the ridges coincide (Method 1/2). Symbol shapes indicate mode degree: I — (circles), / = 1 
(triangles), I = 2 (squares) and I = 3 (diamonds). 



filled symbols show those of a Cen B (Kjeldsen et al. 2005) after multiplying by 1.0478. 
3.1 Measuring mean densities 

Assuming that a given pair of stars are homologous, (i.e., that Equation 1 holds), the value of the scaling 
factor gives a direct measurement of their relative densities, without the need to refer to theoretical 
models. The scaling factors used in Fig. 1 are precise to about 0.05%, in the sense that changing them 
by this amount produces a noticable deviation from parallelism. If we accept the validity of Equation 1 
then our results allow us to measure mean densities for a Cen A and B relative to solar with a precision 
of 0.1%. We obtain (0.8601 ± 0.0003) gcrrr 3 for a Cen A and (2.0018 ± 0.0008) gcrrr 3 for a Cen B. 
These values agree with those found by comparing the observed frequencie of radial modes with models 
that have been corrected for the near-surface offset (Kjeldsen et al. 2008; Table 2), but are more precise 
and do not make any use of model calculations. In practice, however, Equation 1 may not be accurate 
to this level. In particular, defining a mean density requires that we specify the position of the surface, 
which is ambiguous, as discussed by Bahcall & Ulrich (1988) in the context of helioseismology. In any 
case, we have certainly derived mean densities for these stars that are precise enough for any practical 
application. 

The two stars shown in Figure 2, r Cet and a Cen B, are even closer. Once again, we found that 
changing the scaling factor by about 0.05% produced a noticable departure from parallelism. Using the 
mean density found above for a Cen B, the implied mean density for r Cet is 2.198 ± 0.004gcrrr 3 , 
which agrees with the value of 2.21 ± 0.01 gem -3 found by Teixeira et al. (2009) but is more precise. 
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Again, we note that being able to measure the homology scaling factor to high precision does not 
necessarily provide a density measurement with similar accuracy. 

3.2 Ridge identification in F stars 

An important application of scaled echelle diagrams is to the problem of ridge identification in F stars. 
This problem has arisen in the context of several F-type main-sequence stars observed using the CoRoT 
spacecraft. The first and best-studied example is HD 49933 (At/ « 85 /iHz), whose echelle diagram 
from 60 days of CoRoT observations showed two broad and very similar ridges (Appourchaux et al. 
2008). It was clear that one ridge was due to I = and I = 2 modes and the other to I = 1, but the 
combination of significant rotational splitting and large linewidths made it difficult to decide which was 
which. Appourchaux et al. (2008) made a global fit to the line profiles, which led them to favor the 
possibility they labelled 'Scenario A'. Further analysis of the same data has been carried out by several 
groups (Benomar et al. 2009a; Gruberbauer et al. 2009; Roxburgh 2009) and none favored a definite 
identification, while comparison with theoretical models (Kallinger et al. 2010) gave a better match to 
Scenario B. Subsequently, the analysis of an additional 180 days of CoRoT observations using revised 
methods led Benomar et al. (2009b) to reverse the original identification in favor Scenario B. 

Two other F stars observed by CoRoT have presented the same problem, namely HD 181906 
(Af w 87.5 /iHz; Garcia et al. 2009) and HD 181420 (At; w 75 /zHz; Barban et al. 2009). In neither 
case were the authors able to decide the correct scenario. 

Using scaled echelle diagrams, together with the quite reasonable assumption that e varies slowly 
with stellar parameters, we might hope to be able to tie these stars together. Figure 3 shows how this 
works for two CoRoT targets, HD 49933 and HD 181420. In all three panels, the filled symbols show 
Scenario B for HD 49933 (Benomar et al. 2009b). The open symbols show Scenario 1 for HD 181420 
(Barban et al. 2009) with three different scaling factors. In the upper panel, the scaling factor was 
chosen to align the ridges as closely as possible, again using trial and error, and we indeed see a good 
match. However, we should check whether shifting one star by half an order also produces a match. 
This requires changing the scaling factor by 0.5/n max , which is about 2.5% in this case. This is shown 
in the lower two panels of Figure 3, where the scaling factor has been changed in both directions by 
this amount and then fine-tuned to align the ridges. Neither of these match as well as the upper panel, 
giving us confidence that Scenario B for HD 49933 is equivalent to Scenario 1 for HD 181420. That 
is, the two scenarios are either both correct or both wrong. 

The third problematic CoRoT target mentioned above, which has a significantly lower signal-to-noise 
ratio, is HD 181906 (Garcia et al. 2009). The power spectrum of this star is shown as the greyscale 
in Figure 4. Overlaid with filled symbols are the oscillation frequencies for HD 49933 from the revised 
identification (Benomar et al. 2009b; Scenario B) after multiplying by 1.011. There is good agreement 
between the stars and, using HD 49933 as a guide, we are able to follow the I = 1 ridge of HD 181906 
down to quite low frequencies. Examining the two possible identifications proposed for HD 181906 by 
Garcia et al. (2009), we can identify Scenario B for that star with Scenario B for HD 49933. Once 
again, either both are correct or both are wrong. 

Having linked these three F-type CoRoT targets, all of which have quite similar values of At/, 
we would clearly like to confirm the identifications by tying them to other stars whose identifications 
are secure. We do this in Figure 5. The greyscale shows the power spectrum of the CoRoT target 
HD 49385 (Af « 56 /iHz; Deheuvels et al. 2010), for which the I = and 2 ridges are clearly resolved. 
The open symbols show frequencies for r\ Boo (Af « 40 /iHz) after multiplying by 1.390, with mode 
identifications that were verified by three sets of observations (Kjeldsen et al. 1995 2003; Carrier et al. 
2005) Finally, the filled symbols show once again the revised identification for HD 49933 (Benomar 
et al. 2009b; Scenario B), this time multiplied by 0.658. These scaling factors differ from unity by 
more than any others we have considered (since Af covers a bigger range). Despite this, we see good 
alignment of the ridges (Method 1) that gives a consistent picture. Interestingly, HD 49385 shows 
significantly more curvature at the lowest orders than the other two stars. 



6 



Scaled oscillation frequencies and echelle diagrams for comparative asteroseismology 





2400 




- 

■ • 
• • 




1 


A 


-l | - 

A 


27 
26 








□ ■ (9 






A 


25 




2200 




□ ■ «o 
□ to* 






A A 
AA 




24 

23 




2000 




■ o • 






AA 




22 




EC 


^ a 


■o» 




A A 




21 




o 1800 


□ 


• • 




AA 




20 




cr 

£ 1600 


- □ «• 
□ ■ 0» 

□ ■ cm 




AA 
A A 
A A 




19 
18 
17 






□ 


■ • 




A A 




16 




1400 


- ■ 


• 
• 




A A 
A 






15 
14 




1200 




■ • 




A 


A 




13 
12 








I 1 








, 1 


11 























20 




40 


60 


80 










Frequency modulo 85.0 ^Hz 






■ • 


ii 




-, | - 

A 




2400 




1 ' 
■ • 


A 


■ • 






A 


26 






■ • 


A 


• 






A 


35 






• 


A M3 O- 


A ■ « 


□ o 




A 


24 


2200 




■ • A ACQ. 


A * ■ 


□ c 


A 




23 






• ■ A 


A □ O 


A . « 


□ 


A 




22 


2000 


— ■ 


• A 


a a 


A ■ • 





A 




21 


EC 


■ 


«i 


□ AO 


A. « 


□ o * 




20 


o 1800 


■ 


•A 


□ AO 


■A* 
• • 


□ a 
□ a 






19 
18 


^equer 


Aa 
Al 


• 
• 


□ AO 

□ « 


• « 


□ 


AO 




17 


* 1600 


^ A ■ 


• 


□ OA 


■ • 


□a o 




16 




u A ■ 


• 


□ OA 


■ »A 


□ A 


o 






1400 


— ■ 


• 


A 


• 


A 






14 






• 


A 


• 

r 


A 


A 


, 1 


13 
12 


1200 


■ 


• 

1 


A 

A 




) 20 


40 


80 


80 


11 





20 


40 60 80 


Frequency 


modulo 85.0 fjMz 








Frequency modulo 85.0 /zHz 



Figure 3: Scaled echelle diagrams showing two CoRoT F-type targets that have ambiguous identifications. The 
filled symbols are oscillation frequencies for HD 49933 from the revised identification (Benomar et al. 2009b; 
Scenario B). The open symbols are frequencies for HD 181420 (Barban et al. 2009; Scenario 1) after multiplying 
by 1.144 (upper panel), 1.115 (lower left) and 1.173 (lower right). Symbol shapes indicate mode degree: / = 
(circles), I — 1 (triangles), and I = 2 (squares). 
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Figure 4: Scaled echelle diagram comparing HD 49933 with another CoRoT F-type target that has an ambiguous 
identification. The greyscale is the power spectrum of HD 181906 (Garci'a et al. 2009), smoothed to a FWHM 
of 3 ^iHz. The filled symbols are oscillation frequencies for HD 49933 from the revised identification (Benomar 
et al. 2009b; Scenario B) after multiplying by 1.011. Symbol shapes indicate mode degree: I = (circles), 
I — 1 (triangles), and I = 2 (squares). 

To summarise, we conclude that the correct identifications are: Scenario B for HD 49933 (Benomar 
et al. 2009b), Scenario 1 for HD 181420 (Barban et al. 2009) and Scenario B for HD 181906 (Garcia 
et al. 2009). The first two of these agree with the conclusions of Mosser & Appourchaux (2009), which 
were based on autocorrelation analysis of the time series. 

4 Conclusions 

We have described a method for scaling oscillation frequencies and displaying two or more stars on a 
single echelle diagram. Assuming that two stars are sufficiently similar to be homologous, the diagram 
can be used to infer the ratio of their mean densities very precisely, without reference to models. In 
addition, data from the star with the better signal-to-noise ratio can be used to confirm weaker modes 
and reject sidelobes in data from a second star. A very important application is to provide a solution 
to the problem of ridge identification in F-type stars observed by CoRoT, as discussed in Section 3.2, 
and we have successfully applied the method to Procyon (Bedding et al. 2010a). Another application 
is to apply what might be called ensemble asteroseismology to the very large samples of stars being 
observed by the CoRoT and Kepler space missions. The results of applying this technique to red giants 
observed with Kepler are described by Bedding et al. (2010b). 
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Figure 5: Scaled echelle diagram to establish the correct identification for F-type stars. The greyscale is the 
power spectrum of HD 49385 (Deheuvels et al. 2010), smoothed to a FWHM of 1 pHz. The filled symbols are 
oscillation frequencies for HD 49933 from the revised identification (Benomar et al. 2009b; Scenario B) after 
multiplying by 0.658. Open symbols are frequencies for r\ Boo (Kjeldsen et al. 2003) after multiplying by 1.390. 
Symbol shapes indicate mode degree: I = (circles), I = 1 (triangles), and I = 2 (squares). 
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